------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  /Users/caseyklofstad/Dropbox/Effects, Impacts and Influences of Hurricane Irma/irma paper 2 BUSINESS price 
> gouging experiment/replication files/study 1 log.log
  log type:  text
 opened on:   8 Feb 2023, 16:10:29

. 
. *get data
. use "/Users/caseyklofstad/Library/CloudStorage/GoogleDrive-c.klofstad@miami.edu/My Drive/CASEY/_UM_drive/research/_old
> _projects/QAnon MC post/test of irma data.dta"

. 
. *extract vars needed for analysis
. keep Q52a Q52b Q52c Q4 Q3 Q80 Q13 Q14_1 Q14_2 Q14_3 Q14_4 Q14_5 Q14_6 Q10a Q10a

. 
. *create factor variable for experimental condition
. generate gouge=.
(2,085 missing values generated)

.                 *control
. replace gouge=1 if Q52a !=.
(717 real changes made)

. tab Q52a

       Q52a |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        646       90.10       90.10
          2 |         43        6.00       96.09
          3 |         28        3.91      100.00
------------+-----------------------------------
      Total |        717      100.00

. tab gouge

      gouge |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        717      100.00      100.00
------------+-----------------------------------
      Total |        717      100.00

. tab gouge Q52a  

           |               Q52a
     gouge |         1          2          3 |     Total
-----------+---------------------------------+----------
         1 |       646         43         28 |       717 
-----------+---------------------------------+----------
     Total |       646         43         28 |       717 

.                 *stores might run out of goods
. replace gouge=2 if Q52b !=.
(650 real changes made)

. tab Q52b

       Q52b |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        579       89.08       89.08
          2 |         34        5.23       94.31
          3 |         37        5.69      100.00
------------+-----------------------------------
      Total |        650      100.00

. tab gouge

      gouge |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        717       52.45       52.45
          2 |        650       47.55      100.00
------------+-----------------------------------
      Total |      1,367      100.00

. tab gouge Q52b

           |               Q52b
     gouge |         1          2          3 |     Total
-----------+---------------------------------+----------
         2 |       579         34         37 |       650 
-----------+---------------------------------+----------
     Total |       579         34         37 |       650 

.                 *stores close
. replace gouge=3 if Q52c !=.
(717 real changes made)

. tab Q52c

       Q52c |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        621       86.61       86.61
          2 |         36        5.02       91.63
          3 |         60        8.37      100.00
------------+-----------------------------------
      Total |        717      100.00

. tab gouge

      gouge |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        717       34.40       34.40
          2 |        650       31.19       65.60
          3 |        717       34.40      100.00
------------+-----------------------------------
      Total |      2,084      100.00

. tab gouge Q52c

           |               Q52c
     gouge |         1          2          3 |     Total
-----------+---------------------------------+----------
         3 |       621         36         60 |       717 
-----------+---------------------------------+----------
     Total |       621         36         60 |       717 

.                 *summary
. label define gouge_label 1 "control" 2 "stores run out" 3 "stores close"

. label values gouge gouge_label

. tab gouge

         gouge |      Freq.     Percent        Cum.
---------------+-----------------------------------
       control |        717       34.40       34.40
stores run out |        650       31.19       65.60
  stores close |        717       34.40      100.00
---------------+-----------------------------------
         Total |      2,084      100.00

. summarize gouge

    Variable |        Obs        Mean    Std. dev.       Min        Max
-------------+---------------------------------------------------------
       gouge |      2,084           2    .8297169          1          3

. 
. *create create factor variable for experimental condition, this time control vs all treats
. generate gouge_pooled=.
(2,085 missing values generated)

.                 *control
. replace gouge_pooled=1 if Q52a !=.
(717 real changes made)

. tab Q52a

       Q52a |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        646       90.10       90.10
          2 |         43        6.00       96.09
          3 |         28        3.91      100.00
------------+-----------------------------------
      Total |        717      100.00

. tab gouge_pooled

gouge_poole |
          d |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        717      100.00      100.00
------------+-----------------------------------
      Total |        717      100.00

. tab gouge_pooled Q52a   

gouge_pool |               Q52a
        ed |         1          2          3 |     Total
-----------+---------------------------------+----------
         1 |       646         43         28 |       717 
-----------+---------------------------------+----------
     Total |       646         43         28 |       717 

.                 *stores might run out of goods
. replace gouge_pooled=2 if Q52b !=.
(650 real changes made)

. tab Q52b

       Q52b |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        579       89.08       89.08
          2 |         34        5.23       94.31
          3 |         37        5.69      100.00
------------+-----------------------------------
      Total |        650      100.00

. tab gouge_pooled

gouge_poole |
          d |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        717       52.45       52.45
          2 |        650       47.55      100.00
------------+-----------------------------------
      Total |      1,367      100.00

. tab gouge_pooled Q52b

gouge_pool |               Q52b
        ed |         1          2          3 |     Total
-----------+---------------------------------+----------
         2 |       579         34         37 |       650 
-----------+---------------------------------+----------
     Total |       579         34         37 |       650 

.                 *stores close
. replace gouge_pooled=2 if Q52c !=.
(717 real changes made)

. tab Q52c

       Q52c |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        621       86.61       86.61
          2 |         36        5.02       91.63
          3 |         60        8.37      100.00
------------+-----------------------------------
      Total |        717      100.00

. tab gouge_pooled

gouge_poole |
          d |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        717       34.40       34.40
          2 |      1,367       65.60      100.00
------------+-----------------------------------
      Total |      2,084      100.00

. tab gouge_pooled Q52c

gouge_pool |               Q52c
        ed |         1          2          3 |     Total
-----------+---------------------------------+----------
         2 |       621         36         60 |       717 
-----------+---------------------------------+----------
     Total |       621         36         60 |       717 

.                 *summary
. label define gouge_pooled_label 1 "control" 2 "treat pooled"

. label values gouge_pooled gouge_pooled_label

. tab gouge_pooled

gouge_pooled |      Freq.     Percent        Cum.
-------------+-----------------------------------
     control |        717       34.40       34.40
treat pooled |      1,367       65.60      100.00
-------------+-----------------------------------
       Total |      2,084      100.00

. summarize gouge_pooled

    Variable |        Obs        Mean    Std. dev.       Min        Max
-------------+---------------------------------------------------------
gouge_pooled |      2,084     1.65595    .4751715          1          2

. 
. *balance analysis
.         *clean up covariates
.                 *sex
. tab Q4

         Q4 |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |      1,031       49.54       49.54
          2 |      1,050       50.46      100.00
------------+-----------------------------------
      Total |      2,081      100.00

. generate female=.
(2,085 missing values generated)

. replace female=0 if Q4==1
(1,031 real changes made)

. replace female=1 if Q4==2
(1,050 real changes made)

. tab Q4 female

           |        female
        Q4 |         0          1 |     Total
-----------+----------------------+----------
         1 |     1,031          0 |     1,031 
         2 |         0      1,050 |     1,050 
-----------+----------------------+----------
     Total |     1,031      1,050 |     2,081 

.                 *age
. encode Q3, generate(age)

. tab age

        age |      Freq.     Percent        Cum.
------------+-----------------------------------
         18 |         49        2.35        2.35
         19 |         35        1.68        4.03
         20 |         43        2.06        6.09
         21 |         39        1.87        7.96
         22 |         31        1.49        9.45
         23 |         34        1.63       11.08
         24 |         22        1.06       12.13
         25 |         41        1.97       14.10
         26 |         35        1.68       15.78
         27 |         30        1.44       17.22
         28 |         23        1.10       18.32
         29 |         33        1.58       19.90
         30 |         34        1.63       21.53
         31 |         25        1.20       22.73
         32 |         27        1.29       24.03
         33 |         55        2.64       26.67
         34 |         27        1.29       27.96
         35 |         57        2.73       30.70
         36 |         40        1.92       32.61
         37 |         41        1.97       34.58
         38 |         33        1.58       36.16
         39 |         32        1.53       37.70
         40 |         33        1.58       39.28
         41 |         22        1.06       40.34
         42 |         27        1.29       41.63
         43 |         33        1.58       43.21
         44 |         31        1.49       44.70
         45 |         37        1.77       46.47
         46 |         34        1.63       48.11
         47 |         39        1.87       49.98
         48 |         31        1.49       51.46
         49 |         35        1.68       53.14
         50 |         38        1.82       54.96
         51 |         38        1.82       56.79
         52 |         30        1.44       58.23
         53 |         30        1.44       59.66
         54 |         42        2.01       61.68
         55 |         33        1.58       63.26
         56 |         26        1.25       64.51
         57 |         25        1.20       65.71
         58 |         27        1.29       67.00
         59 |         31        1.49       68.49
         5o |          1        0.05       68.54
         60 |         36        1.73       70.26
         61 |         46        2.21       72.47
         62 |         34        1.63       74.10
         63 |         35        1.68       75.78
         64 |         44        2.11       77.89
         65 |         42        2.01       79.90
         66 |         43        2.06       81.97
         67 |         40        1.92       83.88
         68 |         52        2.49       86.38
         69 |         33        1.58       87.96
         70 |         46        2.21       90.17
         71 |         47        2.25       92.42
         72 |         33        1.58       94.00
         73 |         24        1.15       95.16
         74 |         15        0.72       95.88
         75 |         18        0.86       96.74
         76 |         18        0.86       97.60
         77 |         13        0.62       98.23
         78 |          7        0.34       98.56
         79 |          9        0.43       98.99
         80 |          6        0.29       99.28
         81 |          3        0.14       99.42
         82 |          5        0.24       99.66
         83 |          2        0.10       99.76
         84 |          2        0.10       99.86
         85 |          1        0.05       99.90
         88 |          1        0.05       99.95
         89 |          1        0.05      100.00
------------+-----------------------------------
      Total |      2,085      100.00

.                 *income
. tab Q80

        Q80 |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        365       17.51       17.51
          2 |        463       22.21       39.71
          3 |        393       18.85       58.56
          5 |        286       13.72       72.28
          6 |        305       14.63       86.91
          7 |        134        6.43       93.33
          9 |        139        6.67      100.00
------------+-----------------------------------
      Total |      2,085      100.00

. generate income=.
(2,085 missing values generated)

. replace income=1 if Q80==1
(365 real changes made)

. replace income=2 if Q80==2
(463 real changes made)

. replace income=3 if Q80==3
(393 real changes made)

. replace income=4 if Q80==5
(286 real changes made)

. replace income=5 if Q80==6
(305 real changes made)

. replace income=6 if Q80==7
(134 real changes made)

. replace income=7 if Q80==9
(139 real changes made)

. tab Q80 income

           |                                    income
       Q80 |         1          2          3          4          5          6          7 |     Total
-----------+-----------------------------------------------------------------------------+----------
         1 |       365          0          0          0          0          0          0 |       365 
         2 |         0        463          0          0          0          0          0 |       463 
         3 |         0          0        393          0          0          0          0 |       393 
         5 |         0          0          0        286          0          0          0 |       286 
         6 |         0          0          0          0        305          0          0 |       305 
         7 |         0          0          0          0          0        134          0 |       134 
         9 |         0          0          0          0          0          0        139 |       139 
-----------+-----------------------------------------------------------------------------+----------
     Total |       365        463        393        286        305        134        139 |     2,085 

.                 *education
. tab Q13

        Q13 |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |         41        1.97        1.97
          2 |        362       17.38       19.35
          3 |        462       22.18       41.53
          4 |        283       13.59       55.11
          5 |        546       26.21       81.33
          6 |        389       18.67      100.00
------------+-----------------------------------
      Total |      2,083      100.00

.                 *race
. tab Q14_1

      Q14_1 |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |      1,234      100.00      100.00
------------+-----------------------------------
      Total |      1,234      100.00

. tab Q14_2

      Q14_2 |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        404      100.00      100.00
------------+-----------------------------------
      Total |        404      100.00

. tab Q14_3

      Q14_3 |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        552      100.00      100.00
------------+-----------------------------------
      Total |        552      100.00

. tab Q14_4

      Q14_4 |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |         62      100.00      100.00
------------+-----------------------------------
      Total |         62      100.00

. tab Q14_5

      Q14_5 |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |          9      100.00      100.00
------------+-----------------------------------
      Total |          9      100.00

. tab Q14_6

      Q14_6 |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |         10      100.00      100.00
------------+-----------------------------------
      Total |         10      100.00

. generate nonwhite=.
(2,085 missing values generated)

. replace nonwhite=1 if Q14_1==.
(851 real changes made)

. replace nonwhite=0 if Q14_1==1
(1,234 real changes made)

. tab Q14_1 nonwhite

           |  nonwhite
     Q14_1 |         0 |     Total
-----------+-----------+----------
         1 |     1,234 |     1,234 
-----------+-----------+----------
     Total |     1,234 |     1,234 

.                 *party id (based on first screener question)
. tab Q10a

       Q10a |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        829       39.78       39.78
          2 |        589       28.26       68.04
          3 |        539       25.86       93.91
          4 |        127        6.09      100.00
------------+-----------------------------------
      Total |      2,084      100.00

. recode Q10a (1=1) (2=0) (3=0) (4=0), gen(pid_DEM)
(1255 differences between Q10a and pid_DEM)

. recode Q10a (1=0) (2=1) (3=0) (4=0), gen(pid_REP)
(2084 differences between Q10a and pid_REP)

. tab Q10a pid_DEM

           |    RECODE of Q10a
      Q10a |         0          1 |     Total
-----------+----------------------+----------
         1 |         0        829 |       829 
         2 |       589          0 |       589 
         3 |       539          0 |       539 
         4 |       127          0 |       127 
-----------+----------------------+----------
     Total |     1,255        829 |     2,084 

. tab Q10a pid_REP

           |    RECODE of Q10a
      Q10a |         0          1 |     Total
-----------+----------------------+----------
         1 |       829          0 |       829 
         2 |         0        589 |       589 
         3 |       539          0 |       539 
         4 |       127          0 |       127 
-----------+----------------------+----------
     Total |     1,495        589 |     2,084 

.         *multinomial logit with control as base category
. mlogit gouge female age income Q13 nonwhite pid_DEM pid_REP, b(1)

Iteration 0:   log likelihood =  -2281.895  
Iteration 1:   log likelihood = -2271.8218  
Iteration 2:   log likelihood = -2271.8136  
Iteration 3:   log likelihood = -2271.8136  

Multinomial logistic regression                         Number of obs =  2,079
                                                        LR chi2(14)   =  20.16
                                                        Prob > chi2   = 0.1251
Log likelihood = -2271.8136                             Pseudo R2     = 0.0044

--------------------------------------------------------------------------------
         gouge | Coefficient  Std. err.      z    P>|z|     [95% conf. interval]
---------------+----------------------------------------------------------------
control        |  (base outcome)
---------------+----------------------------------------------------------------
stores_run_out |
        female |     .06566   .1173109     0.56   0.576    -.1642651    .2955852
           age |   .0012428   .0034829     0.36   0.721    -.0055836    .0080693
        income |    -.01917   .0360808    -0.53   0.595     -.089887     .051547
           Q13 |  -.0156737   .0433181    -0.36   0.717    -.1005756    .0692282
      nonwhite |  -.0471819   .1321694    -0.36   0.721    -.3062291    .2118652
       pid_DEM |    .045222   .1301807     0.35   0.728    -.2099274    .3003715
       pid_REP |   .0747564   .1431429     0.52   0.601    -.2057985    .3553113
         _cons |  -.0592073    .230816    -0.26   0.798    -.5115984    .3931838
---------------+----------------------------------------------------------------
stores_close   |
        female |   .2383826   .1148428     2.08   0.038     .0132948    .4634704
           age |  -.0031969   .0034152    -0.94   0.349    -.0098905    .0034968
        income |  -.0339373   .0355427    -0.95   0.340    -.1035997     .035725
           Q13 |   .0141893   .0423489     0.34   0.738     -.068813    .0971916
      nonwhite |   .1704618   .1277462     1.33   0.182    -.0799162    .4208398
       pid_DEM |  -.0278272   .1276262    -0.22   0.827    -.2779699    .2223156
       pid_REP |   .1903795   .1399785     1.36   0.174    -.0839732    .4647322
         _cons |  -.0808103   .2259847    -0.36   0.721    -.5237321    .3621115
--------------------------------------------------------------------------------

.                 *note: by adding PID gender now is sig for "stores close"
. 
. *combine and code y/n/dk on PRO anti-gouging law regardless of treatment (include DK responses; 3 = yes ban gouging, 1
>  = no don't, 2 = dk)
. generate no_gouge_all=.
(2,085 missing values generated)

.         *yes, support law
. replace no_gouge_all=3 if Q52a==1 | Q52b==1 | Q52c==1
(1,846 real changes made)

.         *no, do not support law
. replace no_gouge_all=1 if Q52a==2 | Q52b==2 | Q52c==2
(113 real changes made)

.         *dk
. replace no_gouge_all=2 if Q52a==3 | Q52b==3 | Q52c==3
(125 real changes made)

.         *stats/check coding
. tab Q52a

       Q52a |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        646       90.10       90.10
          2 |         43        6.00       96.09
          3 |         28        3.91      100.00
------------+-----------------------------------
      Total |        717      100.00

. tab Q52b

       Q52b |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        579       89.08       89.08
          2 |         34        5.23       94.31
          3 |         37        5.69      100.00
------------+-----------------------------------
      Total |        650      100.00

. tab Q52c

       Q52c |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        621       86.61       86.61
          2 |         36        5.02       91.63
          3 |         60        8.37      100.00
------------+-----------------------------------
      Total |        717      100.00

. tab no_gouge_all

no_gouge_al |
          l |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        113        5.42        5.42
          2 |        125        6.00       11.42
          3 |      1,846       88.58      100.00
------------+-----------------------------------
      Total |      2,084      100.00

. 
. *create indicator whereby 1 = not sure and 0 = yes/no on the ban
. generate dk_gouge=.
(2,085 missing values generated)

.         *yes, support law
. replace dk_gouge=0 if Q52a==1 | Q52b==1 | Q52c==1
(1,846 real changes made)

.         *no, do not support law
. replace dk_gouge=0 if Q52a==2 | Q52b==2 | Q52c==2
(113 real changes made)

.         *dk
. replace dk_gouge=1 if Q52a==3 | Q52b==3 | Q52c==3
(125 real changes made)

.         *stats/check coding
. tab Q52a

       Q52a |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        646       90.10       90.10
          2 |         43        6.00       96.09
          3 |         28        3.91      100.00
------------+-----------------------------------
      Total |        717      100.00

. tab Q52b

       Q52b |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        579       89.08       89.08
          2 |         34        5.23       94.31
          3 |         37        5.69      100.00
------------+-----------------------------------
      Total |        650      100.00

. tab Q52c

       Q52c |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        621       86.61       86.61
          2 |         36        5.02       91.63
          3 |         60        8.37      100.00
------------+-----------------------------------
      Total |        717      100.00

. tab dk_gouge

   dk_gouge |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |      1,959       94.00       94.00
          1 |        125        6.00      100.00
------------+-----------------------------------
      Total |      2,084      100.00

. 
. *combine and code y/n on PRO anti-gouging law regardless of treatment (dropping DK responses; 1 = yes ban gouging, 0 =
>  no don't)
. generate no_gouge_nodk=.
(2,085 missing values generated)

. replace no_gouge_nodk=1 if Q52a==1
(646 real changes made)

. replace no_gouge_nodk=0 if Q52a==2
(43 real changes made)

. replace no_gouge_nodk=1 if Q52b==1
(579 real changes made)

. replace no_gouge_nodk=0 if Q52b==2
(34 real changes made)

. replace no_gouge_nodk=1 if Q52c==1
(621 real changes made)

. replace no_gouge_nodk=0 if Q52c==2
(36 real changes made)

. tab Q52a

       Q52a |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        646       90.10       90.10
          2 |         43        6.00       96.09
          3 |         28        3.91      100.00
------------+-----------------------------------
      Total |        717      100.00

. tab Q52b

       Q52b |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        579       89.08       89.08
          2 |         34        5.23       94.31
          3 |         37        5.69      100.00
------------+-----------------------------------
      Total |        650      100.00

. tab Q52c

       Q52c |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        621       86.61       86.61
          2 |         36        5.02       91.63
          3 |         60        8.37      100.00
------------+-----------------------------------
      Total |        717      100.00

. tab no_gouge_nodk

no_gouge_no |
         dk |      Freq.     Percent        Cum.
------------+-----------------------------------
          0 |        113        5.77        5.77
          1 |      1,846       94.23      100.00
------------+-----------------------------------
      Total |      1,959      100.00

. 
. *ANOVA
. anova no_gouge_nodk gouge female

                         Number of obs =      1,955    R-squared     =  0.0006
                         Root MSE      =    .233539    Adj R-squared = -0.0010

                  Source | Partial SS         df         MS        F    Prob>F
              -----------+----------------------------------------------------
                   Model |  .06002149          3   .02000716      0.37  0.7770
                         |
                   gouge |  .02836244          2   .01418122      0.26  0.7711
                  female |  .03465239          1   .03465239      0.64  0.4255
                         |
                Residual |  106.40852      1,951    .0545405  
              -----------+----------------------------------------------------
                   Total |  106.46854      1,954   .05448748  

. margins gouge

Predictive margins                                       Number of obs = 1,955

Expression: Linear prediction, predict()

---------------------------------------------------------------------------------
                |            Delta-method
                |     Margin   std. err.      t    P>|t|     [95% conf. interval]
----------------+----------------------------------------------------------------
          gouge |
       control  |    .937055   .0089227   105.02   0.000      .919556    .9545539
stores run out  |   .9443527    .009441   100.03   0.000     .9258373    .9628682
  stores close  |   .9455653   .0091224   103.65   0.000     .9276746     .963456
---------------------------------------------------------------------------------

.         *now run interacted w/ pid
. anova no_gouge_nodk gouge##Q10a female

                         Number of obs =      1,955    R-squared     =  0.0065
                         Root MSE      =    .233382    Adj R-squared =  0.0004

                  Source | Partial SS         df         MS        F    Prob>F
              -----------+----------------------------------------------------
                   Model |  .69369213         12   .05780768      1.06  0.3892
                         |
                   gouge |   .0007829          2   .00039145      0.01  0.9928
                    Q10a |  .19964583          3   .06654861      1.22  0.3002
              gouge#Q10a |  .44588676          6   .07431446      1.36  0.2254
                  female |  .01658114          1   .01658114      0.30  0.5812
                         |
                Residual |  105.77485      1,942   .05446697  
              -----------+----------------------------------------------------
                   Total |  106.46854      1,954   .05448748  

. 
. *for supplement: replicate analysis with the "not sure" respondents
.         *not sure vs all other responses
. anova dk_gouge gouge female

                         Number of obs =      2,080    R-squared     =  0.0097
                         Root MSE      =    .236742    Adj R-squared =  0.0082

                  Source | Partial SS         df         MS        F    Prob>F
              -----------+----------------------------------------------------
                   Model |  1.1353154          3   .37843846      6.75  0.0002
                         |
                   gouge |   .6454826          2    .3227413      5.76  0.0032
                  female |  .41897583          1   .41897583      7.48  0.0063
                         |
                Residual |  116.35267      2,076   .05604656  
              -----------+----------------------------------------------------
                   Total |  117.48798      2,079   .05651178  

. margins gouge

Predictive margins                                       Number of obs = 2,080

Expression: Linear prediction, predict()

---------------------------------------------------------------------------------
                |            Delta-method
                |     Margin   std. err.      t    P>|t|     [95% conf. interval]
----------------+----------------------------------------------------------------
          gouge |
       control  |   .0400957   .0088657     4.52   0.000     .0227092    .0574823
stores run out  |   .0574188   .0092941     6.18   0.000      .039192    .0756455
  stores close  |   .0824363    .008853     9.31   0.000     .0650746     .099798
---------------------------------------------------------------------------------

. 
. *save out truncated dataset
. save "/Users/caseyklofstad/Dropbox/Effects, Impacts and Influences of Hurricane Irma/irma paper 2 BUSINESS price gougi
> ng experiment/replication files/study 1 data for replication.dta", replace
file /Users/caseyklofstad/Dropbox/Effects, Impacts and Influences of Hurricane Irma/irma paper 2 BUSINESS price
    gouging experiment/replication files/study 1 data for replication.dta saved

. 
. *close log
. log close
      name:  <unnamed>
       log:  /Users/caseyklofstad/Dropbox/Effects, Impacts and Influences of Hurricane Irma/irma paper 2 BUSINESS price 
> gouging experiment/replication files/study 1 log.log
  log type:  text
 closed on:   8 Feb 2023, 16:10:29
------------------------------------------------------------------------------------------------------------------------
